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ABSTRACT 

Ordinary   differential  equations  with  greatly  different  time 
constants   arise   in  a  wide   variety  of  important  physical  problems.      The 
predictor-corrector  method  of  solution  is   examined  here. 

The   existence   of  stiffly  stable   methods   of  order  as  high  as 
eight  has  been  known  for  some  time.      In  this   report,  we  have   examined 
a  class   of  methods   for  finding  methods   of  order  greater  than  eight. 
We   find  that   the    choice   of 

a(C)   =   5k~r   (S-c)r,    r  =  0,   1,   2,...   k   and  -1  <_  c   <   1 

leads   to  the  methods   of  order  k   as  high   as   eleven   for  some   c   and  r. 
The   multistep  methods   are   stiffly  stable   for  k   <_9   if  r  =   3,    for 
k  <:   10  if  r  =   U,    and  for  k  <_  11  if  r  ■  6. 

We   found  no  twelfth   order  and  higher  stiffly  stable  methods. 
However,   if  we   select 

(Pr)  o(0   =   Ck_r   Ur  -   cr)    or 

k-r     „      _r-i      l 


(q  )      o(?)  =  c~L    z   r  A  c 

i=0 


where   r  =  0,    1,   2,...   k  and  -1  <   c   <   1. 


We   obtain  stiffly  stable  methods   of  P     type   if  k  <_  6   and  of 


Q     type   if  k   <_  7- 
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1 .   INTRODUCTION 

The  mathematical  analysis   of  the  problems   encountered  in 
reactor   calculations,    circuit   analysis   and  chemical  kinetics   often 
leads   to  stiff  systems   of  ordinary  differential  equations,   i.e.    systems 
with  widely  separated  time   constants.      The  standard  methods    for  the 
numerical  integration  of  the  ordinary  differential  equations   are  not  well 
suited  to  the   stiff  differential  equations   for  reasons   of  stability,   the 
time   steps  have  to  be   adopted  to  the   fast   decaying  components  while   one 
is   interested  in  the   slowly   decaying  ones. 

A  number  of  authors   have  proposed  special  algorithms   to 
overcome   this    difficulty   (Lawson    [l],   Richard,   Lanning  and  Torrey    [2], 
and  Liniger   and  Willoughby    [3]).      However,   there  has   remained  still  a 
strong  desire   to   develop  methods   of  Runge-Kutta  or  multistep  type  with 
improved  stability   and  accuracy  properties    so  that  the   above   difficulty 
is    diminished.      The   object   of  this   investigation  is  to   develop  high 
order  stiffly  stable  multistep  methods.      A  linear  k-step  method  for  the 
numerical  solution   of 

y'  =  f(t,  y),  y(tQ)  =  n 

can  be  written  in  the   form 

(1)  k  k 

cu  y     ,   =     £     a     v     ,    .   +  h     E     8.   y'        . 

0  °n+l        .    .      i  °n+l-i  .    .     i  Jn+l-i 

i=l  i=0 

with  a     ^  0  and   |a    |    +    |g    |    ^  0  or  symbolically, 

(2)  e(E)  W*  ■  h  0(E>  ';» 

where  E  is  the  translation  operator   (Ey     =  y        )    and  p    and  a  are 
polynomials   defined  by 

te\  rk  r1^1  rk"2 

p(C)   =   a      C      -   a      C  -   «2   £  -...-   afc 

(3)  a(c)  =  e0  ck  +  e1  c1"1  +  32  ?k"2  +...+  ek 


Therefore,  this  formula  can  only  be  used  if  one  knows  the 
values  of  the  solution  at  k  successive  points.  These  k  values  will 
he  assumed  to  he  given.   Further,  it  can  be  assumed  without  loss  of 
generality  that  the  polynomials  p(0  and  a(£)  have  no  common  factors 
since  in  general  case,  (l)  can  be  reduced  to  an  equation  of  lower  order. 

Definition  1.1.      The  formula   (1)  will  be  said  to  be  of  order  p  >_  0 
if  it  fulfills  the  p  +  1  conditions 

(U)  k 

a  =   I   a. 
°    1=1   X 

k  k 

a  =   Z   (l-i)S  a.  +  s   Z   (l-i)S   6.  ,  s  =  1,  2,...  p 

i=l         x     i=0  X 

Thus  the  method  is  of  order  p  if  for  any  y  e  a  ^         and  for  some 

nonzero  o     ,. 
p+1 

( 5 )  k  k 

a     y         =     I     a.    y     .,    .    +  h     Z      3.    y'.+c,   hP+1  y  PT^ 
0     n+1        .    .      i  Jn+l-i  .    _      i  Jn+l-i  p+1  J(t) 

1=1  1=0 

♦  o(hP+2) 

where  y  "         is   the    (p+l)-st  derivative  of  y  evaluated  for  some  t 
between  t     7  ,    and  t       .      The   last  two  terms  represent  the  truncation 
error. 

Multistep  methods  were   first  investigated  by  Dahlquist    [h] 
who   defined  the   following: 

Definition  1.2.     A  multistep  method  is  called  A-stable  if  all  solutions 
of  (1)   tend  to  zero,   as  n  -»■  °°,   when  it  is  applied  to  the  differential 
equation  of  the  form  y  '  =  Xy  and  X  is  a   (complex)   constant  with  negative 
real  parts. 

He  then  proved  that  the  order  p  of  an  A-stable   linear 
multistep  method  cannot  exceed  two.      The  order  two  method  with  the 
smallest  error  term  is   the  trapezoidal  rule.      Widlund   [5]  has   shown  the 


existence   of  methods   of  orders   three   and  four  which  are   stable   in  the 
wedge    |arg(hX)   —  tt  J    <  a  of  the  negative  half  plane   for  any 

a  <  p-.      Norsett    [6]  has  extended  these  results  to  the  methods   of  orders 

five   and  six.      Gear    [7,8]  has   used  the   conditions  which  are  necessary 
for  stiff  differential  equations.      The   requirements   are   shown  in 
Figure   1. 


STABLE 


hX  =  h  - —  plane 

3y  * 


Figure   1 
Stability  and  Accuracy  Region 


Thus  the  numerical  methods   suitable   for  stiff   differential  equations 
depend  on  the  parameters   D,    0,    a  and  on  the   definition   of  accuracy. 
Gear  has  termed  such  methods   as   stiffly  stable   methods    and  has   obtained 
methods   of  order  as  high   as   six  for  suitable  parameters   D,    0,    and  a. 
Dill    [9]  has   used  computer  graphic  techniques    for  finding  methods   of 
orders   seven   and  eight.      We   shall  show  here   that   the  methods   of  order 
as  high   as   eleven   can  be   obtained  for  suitable  parameters. 


2.      HIGH   ORDER  STIFFLY  STABLE  METHODS 

The   discretization  error  of  the  multistep  method  is   defined 

as  the   difference  between  the  value  y     calculated  from  (l)    and  the 

exact  solution  y(t    ). 
n 

Define 

e     -  y     -  y(t   ) 

n       Jn  n 

Then  the   error  e      obeys   the    difference   equation 


(6) 


(a.    -  hAgJe  =     Z      (a.    +  hAB. )e    .,    .    -  T 

0  0     n+1        .    ,        i  l     n+l-i  n 

1=1 


for  the  equation  y'  =  Ay. 
Assuming  that 

we  get  from  (6)  the  well  known  characteristic  equation  for  £, 

( 7)  k 

(a  -  hX6n)Cn+1  -  £   (a.  +  hAB.)Cn+1-1  +  T  (y,  t  ,  h)  =  0 
0      0        .  ,    i      i  n     n 

i=l 

where  T  is  the  truncation  error  for  the  exact  solution.   For  convergence, 

n 

it  is  necessary  to  bound  the  solution  of  the  imhomogeneous  difference 
equation  (7)  (see,  Henrici  [10]).   It  depends  on  the  stability  of  the 
corresponding  homogeneous  equation  obtained  from  (7)  with  T  set  to  zero. 
This  difference  equation  is  stable  if  and  only  if  all  the  roots  of  the 
polynomial  equation 

k 

(a   -  hAB  )£n+1  -   S   (a.  +  hA6.)C     X  ■  0 
0      0        .  ,    i      i 
1=1 

or  symbolically 

(8)  p(£)  -  hA  o(5)  =  0 


lie   in  or  on  the   unit  circle   and  those   lying  on  the  unit  circle  have 
multiplicity  of  one.      It  should  he  noted  that   it  is   an   asymptotic 
condition  only  and  is   concerned  with  the   convergence   of  y     to  y(t) 
with  h  =    (t   -  t    )/n   as  n  •*■  °°  or  hA  -*■  0.      If  stiff  equations    are  to  he 
integrated  with  large  values   of  h,   then  large  values   of  hX  must  not 
make    (8)    unstable.      Letting  hX  ■>  °°,   the   roots   of   (8)    approach  those   of 
o(C)   =  0.      This   implies   that  the  polynomial  o(^)   must  not  have   roots 
outside   the   unit   circle   and  those   roots    £.    for  which    |£.|    =1  are 
simple.      Further,   a(C)    is    of   degree  k   and  has  no   common  root  with  p(C). 
The   stiff  stability   can  he   investigated  as   follows. 

We  want  hX   values   such  that   (8)   has   roots   inside   the  unit 
circle   or  on  the   unit   circle   and  simple.      The   region   is  hounded  hy  the 
locus   of  p(£)/a(£)    in  the  hX-plane   for   £  =  e      ,    0  e[0s2TT]„      This   locus 
can  be  plotted  in  hX  plane   and  is   of  the   type   as   shown  in  Figure   2. 
At  hX  =   infinity,    the  method  is   stable   if  a(0    is   stable,    so  that   any 
region   connected  to  that  point  will  be   stable  by   continuity   argument.      If 
we  prescribe  k,    then   for  stiffly  stable   methods   the  polynomial  cU) 
should  be   such  that   its   roots    lie  within  the   unit   circle   or  on  the   unit 
circle   and  simple.      Another  restriction   on  aU)    is   that   it  has   no 
common   factor  with  pU).      Since   p(c)   will  always  have   a  root   at  +1, 
a(0   must  not  have   a  root   at  +1.      The   locus    of  p  ( £)  /o(  O    for   C  =  e      , 
0  £[0,2tt]   should  be   as   shown   in  Figure   2. 


Selecting  Coefficients    of  o(E,) 

We   shall   consider  the   stiffly   stable  methods    for  which  k  =  p. 

A  k-step  stiffly  stable   method  of  order  p    (=k)    requires   the  solution 

of  p  +  1  equations   in  2k  +  1  unknowns    (B      can  be  taken  as   l) .      The 

k  +  1  unknowns   a.,   the    coefficients   of  the  polynomial  p(?)    can  be 
l ' 

determined  in  terms   of  k  unknowns    B.  ,   the   coefficients   of  the 
polynomial  c(c)-      The   expressions   for  a.    in  terms   of   B.    as   obtained 
on   solving  the  equations    (k)    are   given  in  Appendix  1. 


Besides  the  above  metnioned  restrictions  on  a(^),  these 

arbitrary  coefficients  8.  can  be  chosen  to 

1 

1.  control  numerical  stability, 

2.  minimize    computational  efforts,    and 

3.  minimize  the  truncation  error. 

Gear  has   chosen   6.    =  0,    i  =   1,   2,...   k,    i.e.    o(£)   =   £      and  obtained 
stiffly  stable   methods   for  k  <_  6.      k  >_  7   does  not   give   a  stable  method. 

These  methods   satisfy    (l)    and   (2).      Dill  has  prescribed  6=1  and 

k— 1 
8     =  -.99,   i.e.    a(0   =   5        (£  -  0.99)    and  obtained  a  seventh  order 

formula.      For  eighth   order  formula  he  has   selected 

_k   2      2 
a(S)   =    £T~    U     -  1.8C  +    .81).      These   formulas   also  satisfy  the 

criterion    (2)  . 

We  have   investigated  a  class   of  methods  which  will  satisfy 
the   essential  property   (l).      Some   of  our  formulas  will  also  satisfy 
(2)    and  (3). 


Class   I 


Choose 


(9)  c(0   =   £k"rU   -   c)r,   r  =  0,    1,   2,...   k 


where 


We   find 


-1   <   c   <   1, 


order  of  stable 
formula 


0 


k<6  k<7  k<8  k<9  k<10  k<ll 


There  are  no  twelfth  order  and  higher  stiffly  stable  methods.  Appendix  2 
contains  some  of  the  special  case  formulas.  The  values  of  the  parameters 
D,    c,   max|^|    (£.    being  the   roots   of  p(C)   =  0)    and  c  (the    coefficients 

of  the   truncation  error)    are   also  given.      The   section  of  the  locus 
p(C)/o(£)    for  seventh  to  eleventh  order  formulas   for   different  values   of 
r  are   shown  in  Figures    3  to   7,    respectively. 


Class   II 


Here  we   take   o(^)    of  the   following  form 


o(0   =   Ck_r  (C 


cr) 


<v 


i=0 
k  and  -1   <   c   <   1. 


whe  re   r  =  0 ,   1 ,   2 , . 

The   above  P      and  Q     type   formulas   for  r  =   0   reduce   to  the 
r  r 

case    considered  by  Gear   and  are   stiffly  stable   for  k   <_  6. 

By  prescribing  the  values   of  c   consistent  with  the   stability 
requirement,    a  locus  was   obtained  which  indicated  the  existence   of 
suitable  parameters   D  and  0  for  a  stiffly  stable  method.      These   results 
are   given  in  the   table  below. 


r 

0              1              2              3              k              5              6 

P 
r 

kl6"     k^T     k^6     k<_7     k<_6     k^6     k<_6 

Q 
r 

k<6     k<7     k<7     k<_6     k<_6     k^6     k<_6 

r  >  6   does  not   indicate   the   existence   of  stiffly   stable   methods. 
Appendix  3  shows  the   coefficients   for  various   order  formulas   of  P 
and  Q     type.      The   values   of  the  parameters    c,   D,   max|c|    and  c  are 

also  tabulated.      The   section   of  the   locus   of  p(£)/o(£)    for  fifth  order 
method  is   shown   in  Figure   8. 


3.   CONCLUSIONS 

The   stiffly  stable  methods   of  order  as  high  as  eight  are 
already  known.      The   aim  of  the  present   investigation  has  been  to   develop 
methods   of  order  higher  than  eight. 

We  have   investigated  the   following  classes   of  methods: 

(1)  oU)  =  Kk"r  (5  -  c)r 

(2)  aU)   =    ?k"r   Ur  -   cr) 

(3)  k-r      r         r    i       i 

0(5)  =  r  r    S    ?       c1 

i=0 

where   r=  0,   1,   2,...    k  and  -1  _f_c   <   1. 

We   find  that   the   choice    (l)    leads  to  the  methods   of  order 
as  high  as   eleven.      The  methods    are   stiffly  stable  for  k   <_  9   if 
r  =    3,    for  k   <_  10   if  r  =   U,    and  for  k   <_  11  if  r  =  6.      Further,  we 
found  no  twelfth   order   and  higher  methods.      The   stiffly  stable   methods 
of  the   types  (2)    and   (3)    do  not   exist   for  k  _>  8.      A  comparitive   study 
of  the  methods    developed  in  this   report   is   under  investigation. 
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APPENDIX  I 


A.   Third  Order  Formula: 


a  v    =   Z  a    v    +  h   £   B  v ' 
0  yn+l   .  .  i+1  °n-i     .  \     i  yn+l-i 
i=0  i=0 


(1.1) 


—        — 

ao 

a 
1 

= 

a2 

a3 

11   2   -1 


18  -3  -6 


-9   6   3  -18 


2-1   2   11 


and 


\   -  -(3S0  -B1+  62-3B3)^  y<*> 


B.      Fourth   Order  Formula: 


0     n+l        .    n     l+l     n-i  .    n     i     n+l-i 

i=0  i=0 


(1.2) 


12 


ao 

al 

a2 

= 

a3 

au 

25 


3-1         1-3 


U8     -10       -8         6     -16 


=     -36       18         0     -18       36 


16       -6         8      10     -U8 


-3         1-1  3       2! 


and 


h5      (5) 


Tn  =  -(12B0-  3B1+2B2-  3B3  t  12  B,,)  gj-^ 


C.      Fifth  Order  Formula: 


0  ''n+l        .    A     l+l  rfn-i  .    ~     l     n+l-i 

i=0  i=0 


(1.3) 


6o 


—       — 

ao 

al 

a2 

= 

a3 

% 

a5 

i 

1 

137       12       -3 


-3         12 


300     -65     -30       15       -20         75 


-300      120      -20      -60  60      -200 


200      -60        60        20     -120        300 


-75       20     -15       30         65     -300 


12 


2        -3         12        137 


—        — 

Bo 

h 

82 

63 

h 

65 

_       _ 

and 


Tn  =  -(10B0  -  2gl  +   S2 


*2tk-mS)hri6] 


D.      Sixth  Order  Formula: 


0     n+l        .    „      i+l     n-1  .    -      i  "n+l-i 

i=0  i=0 


(l.iO 


6o 


*~ " 

ao 

al 

a2 

a3 

= 

% 

a5 

a6 

lUT  10        -2  1-1  2        -10 


360        -77     -2h  9        -8  15        -72 


-U50        150      -35     -U5        30        -50       225 


1+00     -100       80         0     -80       100     -Uoo 


-225  50     -30        U5        35     -150        U50 


72        -15 


_9       2U         77     -360 


-10 


2-1  1-2  10        1U7 


and 


Tn  =  -(60B0   -  103i  +   U32  -   3g3  +   hth  -  1035  +  6066)   |2b-y^T) 


E.      Seventh  Order  Formula: 


7 


0     n+l        .    n     l+l  "'n-i  .    _      l     n+l-i 
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r 

0 

1 

2 

3 

°0 

11 

h9 

257 

15 1+8 

al 

18 

99 

61+2 

1+131 

a2 

-9 

-63 

-531 

-3672 

"3 

2 

13 

lk6 

1089 

60 

6 

30 

200 

2000 

61 

-18 

-280 

-1+200 

62 

98 

29  U0 

63 

-686 

Table   2 
Values   of  the  Parameters 


r 

0 

1 

2 

3 

C 

0 

0.6 

0.7 

0.7 

D 

-0.083 

-0.011+ 

-0.003 

-0.003 

max| c| ,    C¥l 

0.1+26 

0.515 

0.751+ 

0.839 

C.+l 

-0.1361+ 

-0.1837 

-0.3807 

-0.8551 

B.      Fourth  Order  Formulas 


Table   1 
a  -   3   Coefficients 


r 

0 

1 

2 

3 

4 

a 
0 

25 

91 

489 

38849 

5177 

\ 

48 

222 

472 

139934 

1902  4 

a 
2 

-36 

-198 

-1620 

-188922 

-26244 

a 
3 

16 

82 

768 

113330 

l6ll2 

\ 

-3 

-15 

-131 

-25493 

-3715 

60 

12 

48 

300 

32000 

7500 

61 

-36 

-480 

-8l600 

-24000 

62 

192 

69360 

28800 

S3 

-19652 

-15360 

BU 

3072 

Table   2 
Values   of  the  Parameters 


r 

0 

1 

2 

3 

4 

C 

0 

0.75 

0.8 

0.8 

0.8 

D 

-0.667 

-0.164 

-0.018 

-0.007 

-0.011 

max|  £  |  4    £#1 

O.56I 

0.752 

0.821 

0.858 

0.896 

cD+i 

-0.096 

-0.1253 

-0.1849 

-0.35^6 

-0.9740 

Fifth  Order  Formulas 


a  - 


Table   1 
3   Coefficients 


r 

0 

1 

2 

3 

1+ 

5 

ao 

137 

512 

1877 

166528 

22157 

154637 

ai 

300 

1395 

6090 

622875 

92700 

705200 

a2 

-300 

-1560 

-7860 

-925500 

-155460 

-1288700 

a3 

200 

980 

5180 

684500 

130760 

1179 800 

\ 

-75 

-360 

-1815 

-252750 

-55215 

-541175 

a 
5 

12 

57 

282 

37425 

9372 

99152 

6o 

6o 

240 

960 

96000 

15360 

187500 

Bl 

-180 

-l44o 

-216000 

-46080 

-750000 

B2 

5 1+0 

162000 

51840 

1200000 

63 

-40500 

-25920 

-960000 

\ 

486o 

384000 

s 

6l44o 

Table   2 
Values   of  Parameters 


r 

0 

1 

2 

3 

4 

5 

C 

0 

0.75 

0.75 

0.75 

0.75 

0.8 

D 

-2.327 

-0.905 

-0.272 

-O.067 

-0.033 

-0.29 

max|c|  ,    £t*1 

0.709 

0.749 

0.755 

0.793 

0.832 

0.895 

cD+i 

-0.0730 

-0.0898 

-0.1156 

-0.1596 

-0.2507 

-0.7455 

Sixth  Order  Formulas 


a 

Table   1 
-   3  Coefficients 

r 

0 

1 

2 

3 

1+ 

5 

6 

ao 

ll+7 

281+ 

209^ 

7725 

65539 

269927 

0.7519*+ 

a 
1 

360 

797 

7392 

31293 

307832 

13969^0 

1+.  205^8 

a 
2 

-1+50 

-1050 

-11115 

-52965 

-599970 

-3018050 

-9.81212 

"3 

1+00 

900 

9520 

1+861+0 

621520 

31+85600 

12.22517 

ak 

-225 

-250 

-5070 

-26055 

-361265 

-2270525 

-8.57888 

a 
5 

72 

159 

15  81+ 

7875 

111861+ 

-711252 

3.215 

a6 

-10 

-22 

-217 

-1063 

-11+1+1+2 

-115290 

-0.50269 

eo 

60 

120 

960 

381+0 

37500 

187500 

1 

61 

-60 

-lUUO 

-86 1+0 

-120000 

-750000 

-5.1 

B2 

5^0 

61+80 

1U 1+000 

1200000 

10.8375 

33 

-1620 

-76800 

-960000 

-12.2825 

61+ 

15360 

381+000 

7.83009 

65 

-6lUU0 

-2.66223 

^ 

0.37715 

Table  2 
Values   of  the  Parameters 


r 

0 

1 

2 

3 

1+ 

5 

6 

C 

0 

0.5 

0.75 

0.75 

0.8 

0.8 

O.85 

D 

-6.071 

-3.607 

-1.191 

-0 . 1+09 

-0.078 

-0.039 

-0.051 

max| £    ,    5^1 

0.863 

0.789 

0.760 

0.785 

0.837 

0.896 

0.930 

Vi 

-0.0583 

-0.0651+ 

-0.081+3 

-0.1071 

-O.1569 

-0.2521 

-O.8605 
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APPENDIX  III 


A.   Third  Order  Formulas 


Polynomial 


Table   1 

a  -   8 

Coefficients 

1 

P2 

Q2 

P3 

S 

0 

185 

16? 

1321 

1126 

1 

3^2 

29U 

2007 

1962 

2 

-171 

-165 

-639 

-95^ 

3 

lli 

38 

-hi 

118 

0 

96 

96 

750 

750 

1 

0 

-2k 

0 

-U50 

2 

-5U 

6 

0 

270 

3 

-162 

-162 

Table  2 
Values  of  the  Parameters 


Polynomial 

P2 

Q 
H2 

P3 

Q 
3 

C 

-0.75 

-0.25 

0.6 

-0.6 

D 

-0.012 

-0.06 

-O.05U 

-o.oi+o 

max|  £  |  ,    £#1 

0.7^7 

0.1+77 

0.6128 

0.5529 

CP+1 

-0.105U 

-0.1587 

-O.1726 

-0.2558 

B.   Fourth  Order  Formulas 


a  - 


Table  1 
(3  Coefficients 


Polynomial 

P2 

^2 

P3 

Q 
^3 

% 

ao 

101 

93 

1063 

185 

13897 

al 

200 

204 

212.8 

402 

29844 

°2 

-144 

-180 

-188U 

-342 

-23904 

a3 

56 

84 

880 

158 

8812 

% 

-11 

-15 

-61 

-33 

-855 

B0 

48 

48 

500 

96 

7500 

61 

0 

-24 

0 

-48 

-4500 

82 

-12 

12 

0 

2k 

2700 

B3 

256 

-12 

-1620 

"l, 

972 

Table  2 
Values  of  the  Parameters 


Polynomial 

P2 

Q 

^2 

P3 

Q3 

Q4 

C 

-0.5 

-0.5 

-0.8 

-0.5 

-0.6 

D 

-0.557 

-0.4l4 

-0.462 

-0.472 

-0.363 

max] C| ,   C^l 

0.527 

0.662 

0.833 

0.578 

0.6545 

Vi 

-0.0911 

-0.1204 

-0.0820 

-0.1243 

-0.1504 

Fifth  Order  Formulas 


Table 

1 

a   - 

3    Coefficients 

Polynomial 

P2 

^2 

P3 

S 

pu 

% 

s 

ao 

551 

3137 

138U58 

208 

U39 

2077 

U1U2 

ai 

1230 

8320 

310935 

517 

96U 

5150 

10225 

a 
2 

-1180 

-10220 

-3U3740 

-572 

-972 

-5660 

-11120 

a 
3 

7^0 

7160 

21)4580 

388 

66U 

3760 

7220 

\ 

-285 

-2515 

-53130 

-1U8 

-253 

-1U15 

-2530 

a5 

U6 

392 

9813 

23 

36 

2 1+2 

3^7 

eo 

2H0 

1500 

60000 

96 

192 

960 

1920 

61 

0 

-1200 

0 

-1+8 

0 

-U80 

-960 

32 

-6o 

960 

0 

2h 

0 

2 1+0 

1+80 

33 

^37^0 

-12 

0 

-120 

-2U0 

\ 

-12 

60 

120 

S 

-60 

Polynomial 


C 

D 

max|^|  ,    &1 


Table   2 
Values   of  the  Parameters 


Q, 


Q, 


-0.5  -0.8              -0.9  -0.5 

-2.325  -1.077         -1.101  -1.571 

0.637  0.892          0.9^9  0.728 

-0.0708  -0.0975  -0.0670  -0.0875 


-0.5 

-0.5 

-0.5 

-2.250 

-1.553 

_l.i4.U9 

0.7U1 

0.6875 

0.722 

-0.0720 

-0.0886 

-0.0913 

Sixth  Order  Formulas 


a  -    B 


Table   1 
Coefficients 


Polynomial 

Q2 

P3 

s 

?u 

Q 

^5 

ao 

3UU3 

18439 

5655 

1^721+01 

4520 

°i 

10156 

45576 

15655 

3619208 

12493 

a2 

-1U810 

-59130 

-21125 

-4572030 

-16790 

a3 

13280 

50000 

3760 

4192080 

14780 

a4 

-7105 

-252U5 

-10525 

-2334035 

-8200 

a5 

22208 

8424 

3306 

662376 

2615 

a6 

-306 

-1186 

-455 

-95198 

-378 

eo 

1500 

7500 

2400 

6xl05 

1920 

*1 

-1200 

0 

-1200 

0 

-960 

32 

960 

0 

600 

0 

480 

33 

3840 

-300 

0 

-240 

34 

-i44o6o 

120 

35 

-60 

Table   2 
Values   of  the  Parameters 


Polynomial 

Q2 

P3 

Q 
3 

P4 

Q 
^5 

C 

-0.8 

-0.8 

-0.5 

-0.7 

-0.5 

D 

-2.069 

-4.452 

-3.808 

-4.694 

-3.944 

max|  E,\  ,    £^1 

0.934 

0.868 

O.861 

0.933 

0.839 

Vi 

-0.0732 

-O.0566 

-0.0671 

-0.0573 

-0.0677 

E.   Seventh  Order  Formulas 


Table  1 
a  -  3  Coefficients 


Polynomial 

Q 

^2 

P3 

a 
0 

412310 

136U0068U 

a 
1 

13^2600 

370394682 

a2 

-2287481 

-568618092 

a3 

2519300 

605263295 

ak 

-1809150 

-430428180 

a5 

853160 

211815954 

a6 

-234675 

-59320212 

aT 

28556 

7293237 

eo 

168000 

525xl05 

Bl 

-1U2800 

0 

e2 

121380 

0 

33 

28946820 

Table  2 
Values  of  the  Parameters 


Polynomial 

Q2 

P3 

C 

-0.85 

-0.82 

D 

-4.530 

-12.022 

max|c|  ,    £7^1 

0.992 

0.927 

Vi 

-0.0589 

-0.0474 

UO 
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